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We propose a new nonrelativistic Chern-Simons theory based on a simple modification of the standard La- 
grangian. This admits asymptotically nonvanishing field configurations and is applicable to the description of 
systems of repulsive bosons. The new model supports topological vortices and has a self-dual limit, both in the 
pure Chern-Simons and in the mixed Chern-Simons-Maxwell cases. The analysis is based on a new formulation of 
the Chern-Simons theories as constrained Hamiltonian systems. 
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Chern-Simons vortices are currently considered as a 
candidate for anyon-like objects in quasiplanar systems 
of condensed matter physics. Both relativistic and non- 
relativistic matter fields can support these structures [Q . 
So far the relativistic model seems to have exhibited a 
wider spectrum of solutions. It was shown to possess 
both topological vortices, i.e. asymptotically nonvanish- 
ing solutions with quantized flux j2|^] and nontopolog- 
ical solitons for which the matter field ip vanishes 
at infinity. On the contrary, the nonrelativistic model, 
which may happen to be more relevant for actual applica- 
tions, exhibits only asymptotically vanishing solitons, or 
"lumps" j5j. These structures arise when the interaction 
between the bosons is attractive, U ~ —(ipip*) 2 . Since 
the lumps have nonzero vorticity and even some topo- 
logical characteristics (the flux is quantized), they are 
sometimes referred to as vortices and topological solitons. 
However, these are clearly distinct from what are usually 
called topological solitons on the plane, i.e. asymptoti- 
cally nonvanishing solutions approaching different vacua 
in different directions. The "genuine" topological vor- 
tices, or "topological defects" interpolating between dif- 
ferent nonzero vacua, turned out to be absent in the non- 
relativistic Chern-Simons theory, though one could have 
expected them to arise in the case of repulsive nonlinear- 
ities U(ip). Furthermore, it can be readily shown that 
the nonzero vacuum solution itself is not admitted by the 
"standard" nonrelativistic model. There is no "conden- 
sate" in this theory and it cannot be used to describe 
repulsive boson gases. 

The main objective of the present work is to incorpo- 
rate the condensate into the nonrelativistic gauge the- 
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ory. We show that the inapplicability of the standard 
model for the description of asymptotically nonvanishing 
fields can be traced back to an inadequate Lagrangian 
formulation of its nongauged precursor. Appropriately 
modifying this Lagrangian and using it as a basis for the 
gauge theory, we arrive at a new version of the gauged 
nonlinear Schrodinger equation which is completely com- 
patible with the nonvanishing ( "condensate" ) boundary 
conditions. 

We start with the one-dimensional nonlinear Schro- 
dinger equation with a general nonlinearity, 

irh+^xx + F(fi)rl> = 0. (1) 

Here p = |?/>| 2 , and F(p) = —dll/dp. The standard La- 
grangian for eq. (|l|), 

c= l -{vj t r-r t i>)-\^\ 2 ~u( P ), (2) 

does not automatically produce correct integrals of mo- 
tion for solutions with \%j]\ 2 approaching po at infinity. 
First of all, the integral for the number of particles, 

takes the form N = J p dx and apparently diverges. The 
regularized number of particles, 

N = J (p-po)di, (4) 

is obtained by the ad hoc subtraction of the background 
contribution. 

The standard definition of momentum, 

does not yield the correct expression either. For the La- 
grangian (^) one obtains P = (i/2) J (ip x ?p* — xjj*xjj)dx 
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which is not compatible with the Hamiltonian structure 
of the model M. Indeed, varying this P gives 



5P = i (ipxSip* ~ ipx5ijS)dx ~ PoSArgip 



+00 



Due to the appearance of the boundary term we cannot 
compute the functional derivatives SP/Sip and SP/Sip*. 
As a result, the Poisson bracket of P with some other 
functional (e.g. the Hamiltonian) can not be evaluated. 
To obtain the definition compatible with the Hamiltonian 
structure we again have to make an a posteriori regular- 
ization @] (see also ||): 



P 



ipip* x ) dx + p Argip 
Po 



tyv* - Wl) ( 1 - 



+00 



p 



dx. 



(6) 



Proceeding to two dimensions, the standard definition 
P = (i/2) J (ip* Wip ^ ipVip*) d 2 r is even less suitable, 
since for asymptotically nonvanishing configurations this 
integral, in general, diverges. 

Now we turn to the gauged nonlinear Schrodinger 
equation as formulated by Jackiw and Pi Q : 



iipt - eA iP + DV + F(p)ip = 0. 



(7a) 



Here D a 
satisfies 



d a + ieA a and the Abelian gauge field A c 



pd F? c 



K 

r 



01 



eJ° 



(7b) 



Eq. (7b) is the most general linear gauge-field equation. 
It combines both the Maxwell and Chern-Simons interac- 
tions. J a — (Jo, J) denotes the conserved matter current: 



J = 7 {V>*(D^)-V>(DV)*}. 



(8a) 



(8b) 



Fa/3 = d a Afj — d a Ap and Greek and Latin indices run 
over 0,1,2 and 1,2, respectively. The parameters p and 
k control the relative contributions of the Maxwell and 
Chern-Simons terms in the corresponding Lagrangian: 



C = - {V* (Dot) - ip(D ip)*} - (D k iP)*(D k TP) 

~^F aP F^ + ^A,F a0 -U( P ). (9) 

p is > 0; we also assume k > 0. The case of negative 
k is recovered by the parity transformation (x 1 ^ x 2 , 
A 1 ^ A 2 ). 

This system inherits the drawback of its nongauged 
counterpart. The momentum, defined by 



Pi = 
equals 

Pi = 



dip dC + dip* dC , dA a dC 



dx 1 dipt dx 1 dip* dx 1 dAf 



d 2 x, (10) 



d 2 x 



= J p{eA l - d t krgiP)d 2 x- p J e l:j E j Bd 2 x. (11) 

Here E 1 — F M and B = F 21 are electric and magnetic 
fields, respectively. The momentum in this case is finite 
provided A 1 — » (l/e)(9iArg ip at infinity, but the func- 
tional derivatives SPi/Sip and SPi/5ip* still cannot be 
calculated. 

Returning to the one-dimensional nongauged nls (Q), 
a natural question is: Is it possible to find a Lagrangian 
producing the same Euler-Lagrange equations as eq. (Q) 
but at the same time automatically yielding the correct 
integrals of motion? The answer is yes. The Lagrangian 



c = Uip t r-r t ^) (1-- 

2 v p 



U(p), (12) 



produces the same nls ([j]), while formulas (||) and (|^) 
yield the regularized integrals N and P, as given exactly 
by (@) and @. 

One may argue that the modification is in a sense triv- 
ial, since eqs. (|l2]) and (^) differ only by a total derivative. 
However, this is no longer the case if the gauge field is 
added: 

C = l - {V*(A)V0 - i>(D m (i - - (Dkip)*{D k ip) 
-^F aP F^ + ^A,F a p-U( P ). (13) 
This is a new system. The field equations resulting 



from (13) differ from those of the standard model (|9j). 
The difference lies in the definition of the number den- 
sity: 



Jo= P- Po- 



(14) 



Apart from that, the field equations are the same, 
eqs. (0), with the spatial part of the current being given 
by the standard expression (8b). However, the new sys- 
tem does possess the condensate solution, while the La- 
grangian (|l^) produces the correct number of particles, 
N = f(p — Po) d 2 r, and the momentum which is compat- 
ible with the Hamiltonian structure of the model: 



P 



l{^(AV)-V>(A#*}(i-- 
2 VP 



-pe^E'B 



d 2 x 



= J[(P- Po) i eAi - 9iAvgip) - lUit&B] d 2 x. (15) 
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Physically, eq. (O) implies that the condensate is elec- 
trically neutral. This may result for instance from the 
presence of motionless particles of the opposite charge. 

The Hamiltonian formulation of this model deserves a 
special comment. For the nonvanishing Maxwell term 
(fi 7^ 0), it can be formulated only as a constrained 
Hamiltonian system. The Hamiltonian reads 



H = 



\^\ 2 + U(p) + ^^ 2 + B 2 ) 



d 2 x 



(16) 



and the constraint is the time component of eq. (7b) (the 
modified Gauss law): 



p divE — kB — e(p — po). 



(17) 



The first pair of canonically conjugate variables is tf> and 
-0*, as usual. The other two canonical coordinates are 
A 1 and A 2 , with the corresponding conjugate momenta 
being 



n, 



dC 



d(d A l ) 2 lJ 



= ZeyA* - pF« 



Eqs. (Q) can be represented then as 
where H is the Lagrange function 



SH 
6IL 



n, 



(18) 



SH 



H = H + A {e(p - po) - (i divE + kB} d 2 x, (19) 



and A — A (x,t) plays the role of the Lagrange multi- 
plier. This formulation remains valid, of course, for the 
standard model (po — 0.) In the pure Chern-Simons case 
(fx = 0) the constraint (|lj]) can be explicitly resolved 
which produces an alternative, unconstrained formula- 
tion g. 

The condensate solution has the form: ip = i/p , 
Ao = (l/e)F(p ), A = 0. This solution exists even when 
the potential U(p) does not possess a symmetry breaking 
minimum at po- We can however confine ourselves to po- 
tentials with U'(po) = as this condition can always be 
accomplished by the transformation Ao — > Ao + F(po)/e, 
U(p) — > U(p) — F(po)p. The condensate solution is then 
ip = y/p , A = 0, A = 0, and it corresponds to an 
extremum of U(p). Through the gauge invariance, the 
condensate generates a set of singular solutions, 



i> = Vf>o e 



iill) 



Ao = 0, A = 



n eg 



(20) 



which serve as asymptotes for vortices as r — > oo. The 
magnetic flux of the vortex is quantized: j B d 2 x = 
<f>Agrd0 = (2n/e)n. Integrating eq. (|l7|) over the en- 
tire plane we observe that the flux and the charge of the 
vortex Q — e J(p — po) d 2 x are related: — k J B d 2 x = Q. 
Hence, Q is quantized as well: 



Q = e J (p - po) d 2 x = — ^^n. 



(21) 



Confining consideration to radially symmetric config- 
urations, we write 



ip(x) = iP(r)e ine 
A (x) =A (r), 



A J (x) 



e jk x k 



$(r) 



3,k = 1,2. 



(22a) 
(22b) 

(22c) 



The magnetic field B(r) — (l/r)<E> r . Regularity at the 
origin requires ip(0) = when n / 0. For finite A(0) 
we should also have $(0) = 0. Now making the gauge 
transformation Arg?/> — > Axgip — n6, $ — > $ — n/e, the 
system (^) is cast in the form 

A?/> - e 2 —^ - eAoi> + F(ip 2 )ijj = 0, 
dA 

pB r 



(23a) 



- 2e 2 — $ = 0, 



dr r 

1.2 



pAA + kB + e(ip z - po) = 0. 



(23b) 
(23c) 



This system does not contain n explicitly. Vortices with 
different vorticities will be distinguished by the boundary 
condition at the origin: $(0) = —n/e. Eq. ( pl| ) indicates 
that conventionally shaped topological vortices for which 
p(r) grows monotonically from to po, are of positive vor- 
ticity. Vorticity-free solutions and those with n < have 
to be nodal, i.e. p(r) — p should necessarily change sign. 
Relegating the comprehensive discussion of solutions to 
the system (^) to a more detailed publication |1(J , here 
we consider only an important special case. 

Static solutions correspond to stationary points of the 
Hamiltonian ([l6|) on the constraint manifold (|l7|) . In the 
mixed CS-Maxwell case, using the flux-vorticity relation, 



jBd 2 



(2n/e)n and the Bogomol'nyi decomposition, 



|D</f 



(£>! ± iD 2 )ip\ 2 i-VxJi eBp, 



(24) 



the energy ( |l6| ) takes the form 



H 




|0D 1 ±iD 2 )V| 2 + ^ 



U[p) 



B±-(p 
A 6 



Po) 



(VA ) 



± -V x J 

2 



d 2 x 
(25) 



For U(p) — (e 2 /2p)(p — po) 2 (which corresponds to the 
Bose gas with the (^-function pairwise repulsion) and 
fields approaching the condensate background (^0|), the 
energy can be rewritten as 



H 



\(D 1 ±iD 2 )ip\ 2 + % 
'(V A ) 2 } d 2 x ±2wp n. 



B±-(p 



Po) 



(26) 
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The lower bound of energy, H = ±2wpon, is saturated 
when the following self-duality equations are satisfied: 



(£>i ± 1D 2 )V> = 0, 

B±Up-p )=< 
M 

A =0. 



(27a) 
(27b) 

(27c) 



The upper (lower) sign should be associated with the 
positive (negative) vorticity n. This can be concluded, for 
example, simply from the fact that the energy is positive 
for positive U(p), see eq. (|l6j). Comparing (|27b|) to (17), 
we see that solutions of eqs. (|^) lie on the constraint 
manifold only for k = p and only in the case of the 
upper sign. Thus, for k > 0, only vortices with positive 
vorticity may exist. Eq. (|27a|) yields 



1 



1 



±—eijdj lnp + -diAvg^p, 



(28) 



and we should retain only the upper sign here. Substi- 
tuting this into (27b), we arrive at 



V 2 lnp = 2-(p- p ) 

K 



(29) 



Eq. ( p9|) appeared previously in the self-dual limit of 
the relativistic Higgs model with Maxwell term and 
is known to possess solutions with "topological vortex" 
asymptotic behaviour: p(oo) — po, p(0) ~ r 2n , n > 1. 
The phase space analysis shows that it also admits lump 
solitons vanishing at infinity [fiof , but these are extrane- 
ous to the subject of the present work. 

The self-duality reduction is also possible in the pure 
Chern-Simons case (p = 0). Making use of the iden- 
tity ( pi| ) and the constraint (|l7j), the energy (16) can be 
represented as 



H 



(A ± iD 2 ) 



T —(p-po)p + U(p) 



Invoking eq. (|21|), we rewrite this as 



H = 



KD^iD^T-ip-poV + Uip) 



d 2 x. 



d 2 x 



±2irpQn. 



(30) 



With the choice of U(p) — ±(e 2 / K)(p — p ) 2 one observes 
that the energy is minimal provided (D\ ± iDijifj — 0, 
whence we have, as before, eq. (p8|). Substituting this 
into the constraint equation (|l7j), we arrive at 



V 2 lnp=±2— (p-p ). 

K 



(31) 



Note that no equation for the scalar potential Aq arises 
here. This is not surprising as Aq is not a dynamical 
variable. (In our Hamiltonian formulation, it is just a La- 
grange multiplier.) For any combination of the Maxwell 



and Chern-Simons terms Aq can be determined from the 
spatial part of eq. (fFi]) |ll[] : 

pV x B + kV x A = eJ. (32) 

The matter current is calculated from eq. 
p. Substituting this into eq. (32) for p = 
solve for Aq: 

Po)- 



Aq = T-(P" 

K 



(|2J): J = T Vx 

0, we can readily 

(33) 



In contrast to the mixed CS-Maxwell case, both signs 
are allowed in ( pl| ) and so the pure Chern-Simons model 
exhibits a wider variety of self-dual solutions jl(| . Here 
we only dwell on the repulsive case, U = (e 2 / ' n){p — po) 2 
(positive sign in (|3l|)). Our remark below pertains to the 
structure of the nonrelativistic Chern- Simons vortices. 
As we have already mentioned, the "vortex" solutions of 
eq. (^9|) are conventionally shaped: p(r) is ~ r 2n near 
the origin, and exponentially approaches p as r — > oo. 
Eq. ( p"7j ) with p = implies that the magnetic field is 
greatest at the centre of the vortex, and rapidly decays 
to zero as r— >oo. (This is in contrast to the relativistic 
case [||J|] where the magnetic field is concentrated in a 
ring.) Eq. ( |33| ) meanwhile indicates that the electric field, 
E = — e r dAo/dr, is localized within a ring around the 
centre of the vortex and vanishes both at r — and r — * oo. 
The structure of the mixed CS-Maxwell self-dual vortex 
is the same, with the exception that it carries no electric 
field. 
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